Collective Synchronization Induced by Epidemic Dynamics 
on Complex Networks with Communities 
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Much recent empirical evidence shows that community structure is ubiquitous in the real-world 
networks. In this Letter, we propose a growth model to create scale-free networks with the tunable 
strength (noted by Q) of community structure and investigate the influence of community strength 
upon the collective synchronization induced by SIRS epidemiological process. Global and local 
synchronizability of the system is studied by means of an order parameter and the relevant finite- 
size scaling analysis is provided. The numerical results show that, a phase transition occurs at Q c ~ 
0.835 from global synchronization to desynchronization and the local synchronization is weakened 
in a range of intermediately large Q. Moreover, we study the impact of mean degree (k) upon 
synchronization on scale-free networks. 
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/. Introduction. — The study of networked systems, in- 
cluding technological, social and biological networks of 
various kinds, has attracted much attention in physics 
community p], [E IE 0j- How the properties of networks, 
such as the lengths of shortest paths between vertices, 
degree distribution, clustering coefficient, degree-degree 
correlation and so on, affect dynamicalprocesses taking 
place upon the networks [E |g, 0, IE IE EE]i nas been 
one of the most important subjects of the body of work. 
Recently, it has been determined that many real-world 
networks show community structure (ill. [T3|. i.e., groups 
of vertices that have a high density of edges within them, 
while a lower density of edges between groups. However, 
there's few work about the influences of various degree 
of community structure upon dynamics. 

In this paper, we intend to fill this gap by investigating 
synchronization behavior induced by the SIRS epidemio- 
logical dynamics [l3j, LLJ] on the scale-free networks with 
various strength (noted by Q) of community structure. 
In Ref. [l5| , the authors have studied the SIRS on small- 
world networks and found that when p, which character- 
izes the degree of disorder of the network, reach an inter- 
mediately large value p c , synchronization of the system 
emerges. Comparatively, we focus on global and local 
(inside each community) dynamics, and discover that no 
synchronization comes forth when the network possesses 
strong enough community structure, i.e. the communi- 
ties are connected by few edges among them. Moreover, 
the vertices inside each community behave weaker syn- 
chronization when Q is in a range of intermediately large 
values. 

II. Network Model. — To generally characterize the 
community structure of scale- free networks, we propose a 



growth model to create a network with a tunable param- 
eter denoting the strength of community structure. In- 
spired by two ingredients of Barabasi- Albert model (BA 
for short), i.e., growth and preferential attachment [16j . 
the rules of our model are as follows: Starting with c com- 
munities, noted by U\, U2, U c —i, U c , and each commu- 
nity with a small number (mo) of vertices. At every time 
step, we add into each community a new vertex with 
m(< mo) edges that link the new vertex to n different 
vertices in this community and m — n different vertices 
in other c — 1 communities already existed in the sys- 
tem. The initial mo x c vertices link to each other to 
keep the connectivity of the network. The values of m 
and n are not necessary integers (take m for example: 
the fractional part of m denotes the probability to link 
m! + 1 different vertices, where ml is the integral part 
of m). When adding a new vertex into community Ui, 
firstly choose n different vertices in community Ui ac- 
cording to "preferential attachment" , which means the 
probability Y[ that the new vertex will connect to ver- 
tex i (i G U{) depends on the degree fcj of vertex i, such 
that n(^i) = kil^Zjeu kj- Then for each one of the 
other m — n edges of the new vertex, choose a commu- 
nity Uh{j^ Ui) randomly and connect the new vertex to 
one vertex in Uh following the preferential attachment 
mechanism referred above. 



The scaling behavior of the degree distribution can be 
calculated by using several approaches [13, EE EH- In 
our model, the degree distributions p(k) of vertices of 
the global network, as well as the local vertices (inside 
each community), are power-law with exponent 3.0, i.e., 
p(k) oc A; -3 (see Fig.l). The analytic procedure is sim- 
ple and not shown here. 
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As proposed by Newman and Girvan [20] and modified 
by Kashtan and Alon [2l|, the strength of community 
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FIG. 1: The global (left) and local (right) degree distribution 
of the network with N = 10 5 , c = 10, m = 4.0 and n = 3.0, 
that says Q = 0.65. It is worthwhile to point out that, for 
different values of Q, the distributions do not change. 



structure can be quantified by 
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where c is the number of communities, L is the number of 
edges in the network, l a is the number of edges between 
nodes in community U s , and d s is the sum of the degrees 
of the nodes in community U s . Roughly speaking, Q is 
the ratio of the number of edges intra-community to the 
total number of the edges. Obviously, if the network is 
divided into some communities more clearly, i.e. there 
arc fewer edges among different communities, the value 
of Q is larger. In our model, for large N (the number 
of all vertices), L — mN, l s = — and d s = 2n + (m — 
n) + (m-n)*(c-l)/(c-l) = Substituting these 

results into Eq. {!]), we obtain 



n 
rn 



(2) 



Thus, for fixed m and c, we modulate the value of n to 
get the networks with various community strength Q. 

III. Epidemic Model. — We analyze SIRS epidemic 
model and aim to point out the role of community struc- 
ture on the temporal dynamics of the epidemic spreading. 
The disease has three stages: susceptible (S), infected (I), 
and refractory (R). A vertex of the networked population 
is described by a single dynamical variable adopting one 
of these three values. Susceptible elements can pass to 
the infected state through contagion by an infected one. 
Infected elements pass to the refractory state after an 
infection time T% . Refractory elements return to the 
susceptible state after a recovery time Tr. The conta- 
gion is possible only during the S phase, and only by an 
I element. During the R phase, the elements are immune 
and do not infect. The system evolves with discrete time 
steps. Each vertex in the network is characterized by a 
time counter Tj(t) = 0, 1, Tj + Tr = T, describing its 
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FIG. 2: The time series of the fraction of infected vertices. 
The systems have N = 10 4 ,c = 25 and m — 4.0, and the 
infection cycle with Ti = 8 and Tr — 5. The left three figures 
(al), (bl) and (cl) showed the global fluctuations of 7ii n f(t) on 
the network with Q = 0.46, 0.81 and 0.935 respectively. The 
right three figures showed the local fluctuations correspond- 
ingly. It's obvious that the global and local fluctuations are 
very different. The detailed analysis is presented in the text. 
The two bottom figures show the clear global (dl) and local 
(d2) periodic oscillations on the network with weak commu- 
nity structure (Q = 0.46). The time steps have been scaled 
by the natural period T of the infection cycle. To is the pe- 
riod of the oscillations. It is manifest that To > T, which is 
different from the result To = T presented for SW networks 
in Ref.fH. 



phase in the cycle of the disease. The epidemiological 
state TTi (S, I, or R) of the vertex depends on the phase 
in the following way: 



Tn(t) = S ifn(t) = 

TTi(t)=I if Ti(t) € [l,Tj] 

TTi(t) = R if Ti(t) 6 + T] 



(3) 



The state of a vertex in the next step depends on its 
current phase in the cycle, and the state of its neighbors 
in the network. A susceptible vertex stays as such, at 
t = 0, until it becomes infected. Once infected, it goes 
(deterministically) over a cycle that lasts T time steps. 
During the first Tj time steps, it is infected and can po- 
tentially transmit the disease to a susceptible neighbor. 
During the last Tr time steps of the cycle, it remains in 
state R, immune and not contagious. After the cycle is 
complete, it returns to the susceptible state. As men- 
tioned in Ref . [l5| , if vertex i is susceptible and it has ki 
neighbors, of which fcj n f are infected, then, i will become 
infected with probability k- m {/ki. 

IV. Results and Analysis. — Specifically we study the 
behavior of the infected sites with respect to Q. A typi- 
cal realization starts with the generation of the network 
characterized by Q and the initialization of the states 
of the vertices. The initial fraction of infected vertices 
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FIG. 3: (Color online) (a) shows the order parameter for 
global synchronization plotted as oN^I" with (3/u = 0.25 vs 
Q for different network size N with fixed Nc = 200, where Nc 
is the number of vertices in each community. There is given a 
unique crossing point at Q c — 0.83(5). From (b) we obtained 
(1 - /3)/V = 0.57(2). These yield /3 » 0.30 and V » 1.22. (c) 
displays the order parameter for local synchronization vs Q 
for different network mean degree (k) — 6, 8, 10 (from bottom 
to top). 




FIG. 4: (Color online) (a) shows the order parameter vs m 
for different natural period T = 8, 10, 12, where m = (k)/2. 
(b) displays the order parameter a for natural period T = 10 
plotted as aN 13 ^ with (3/v = 0.16 vs m for different network 
size N. There is given a unique crossing point at m c = 2.80. 
From (c) we obtained (1 - 0)/V = 0.35(0). These yield w 
0.31 and V « 1.95. 



^inf(O) = 0.1 and the rest susceptible, was used in all the 
simulations here. 

After a transient period, a stationary state is achieved. 
We find that the pronounced fluctuations of the fraction 
of infected vertices is a function of time. Figure 2 shows 
three time series displaying the fraction of infected ver- 
tices in the network with varying community strength 
Q. When Q = 0.46 (see Fig.2(al)), the network has a 
weak strength of community structure. It is similar to 
the real-world networks where the community strength 
Q falls in the range from about 0.3 to 0.7 [2fj]. In such 
condition, the fraction of infected vertices exhibits large 
amplitude oscillations. For strong community structure, 
such as Q = 0.81,0.935 (see Fig.2(bl) and (cl) respec- 
tively), the time series have regular periods but the am- 
plitudes is small and disordered. In addition, we study 
the local dynamics, that is the epidemic process inside 
each community. Figure 2(a2), (b2) and (c2) show the 
time evolution of the fraction of infected vertices in a 
community, for Q = 0.46,0.81 and 0.935 respectively. 
The amplitudes are almost the same. Since the ampli- 
tude is related to the synchronization of the system, we 
will give a measure below and make it clearer. Fig. 2(dl) 
and (d2) show the clear periodic oscillations of the frac- 
tion of local and global infected vertices vs the scaled 
time t/T while Q = 0.46, respectively. One can see that 
the period To is larger than the natural period T of the 
infection cycle, which is different from the result on small- 
world networks presented in Ref. 22]. Moreover, we have 
done the Fourier power analysis for different Q and find 



that there is a sharp peak in the frequency 1/Tq (not 
showed here), which reveals that the series have regular 
temporal periods although the amplitudes are variable. 
As the difference between To and T is the time staying at 
the state S and T is of the same value for all vertices, we 
could analyze the reason of regular To by estimating the 
probability Pi (t) of a vertex i changing state from S to I 
at time t. We here obtain Pi(t) by using the mean-field 
estimation as the following, Pi(t) = ki n f/ki oc rii n f(t). 
That results implies that all the vertices update their 
states at almost the same time which induces the regu- 
lar temporal cycles. Besides the parameter-free infection 
mechanism, there may be other reasonable choices. For 
example, if the susceptible had a probability A of con- 
tagion with each infected neighbor, then the probability 
of infection is [1 — (1 — A) fein/ ]. For small A, we have 
Pi(t) = [1 — (1 - A) fc "^] w A * k in f cc n m f(t). Obviously, 
that dose not affect the qualitative results here. 

To quantify the amplitudes of the oscillation series, we 
define the relevant order parameter 



a(t) 



1 



(4) 



where <fij = lixijj — 1)/T is a geometrical phase corre- 
sponding to Tj. The states r = have been left out of 
the sum in Eq.((4]). We obtain the synchronization order 
parameter a by averaging over 10 4 time steps after the 
transient to a stationary state and subsequently by aver- 
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aging over 400 different realizations of the system. Here 
the synchronization is a measure of the collective order. 
If at any time t all the vertices are almost at the same 
state, i.e. T\ is equal to the same value for all i, the system 
is synchronous. While the vertices are at different states 
equally, the system is not synchronous. Obviously, when 
the system is not synchronized, the phases are widely 
spread in the cycle and the complex numbers e z< ^ are cor- 
respondingly spread in the unit circle which leads to low 
value of a. In contrast, when a significant part of the ver- 
tices are synchronized in the cycle, a is large. The full 
synchronization, i.e. a — 1, will be achieved only when 
all the vertices enter the same state simultaneously. For 
the local synchronization, we calculate the above order 
parameter over the vertices in one community. 

Precise calculation of the critical community strength 
Q c separating synchronized and desynchronized states re- 
quires considering the finite-size effect. In the thermo- 
dynamic limit, the order parameter displays the critical 
behavior a ~ (Q — Q c )^ , with the critical exponent B. 
While in a finite system with size much larger than the 
additional length scale, the critical scaling from is 

<j = N-P' v F[{Q-Q c )N 1 / v l (5) 

where the exponent v describes the divergence of correla- 
tion volume at Q c , ~ \Q — Q c \~ v ■ Since at Q = Q c 
the function F in Eq. © has a value independent of N, 
plotting aN^l v vs Q for various sizes, one can get the 
value of P/V that gives a unique crossing point at Q c . 
One then use 

, r da , 1 — B , „ T 

ln[— ]q c = — In N + const (6) 

in order to determine the value of (1 — B)/u. Then the 
exponent 8 and 77 can be figured out. Figure 3(a) displays 
the determination of Q c for the global synchronization 
using the finite-scale form Eq. (5). Varying the value of 
B/v we find that B/v ~ 0.25 gives a well-defined crossing 
point at Q c ~ 0.835. In Fig. 3(b), the least-square fit 
to Eq. (6) gives (1 - 0)/V ~ 0.572. These yield 8 w 
0.30 and V w 1.22. Moreover, as showed in Fig. 3(c), 
for different mean degree (k), the local synchronization 
parameter falls into the pit around a value Q ~ 0.75. 
This implicates that, when the communities are almost 
unattached (i.e. for very large Q) the local dynamic lies 
on the inner structure of community independently, while 
the communities couple each other strongly (for small Q), 



the local dynamic is almost the same as the global one, 
and in the midst the local synchronization is the weakest. 

Further more, we have studied the impact of the mean 
degree (k) of scale-free networks upon the synchroniza- 
tion. We start the simulation with a generation of scale- 
free networks with Q = 0, that is the BA model or our 
model with n/m = 1/c, and the initial fraction of in- 
fected vertices riinf(O) = 0.1. We let the number of ver- 
tices m that a new added vertex will connect be real 
number, as referred in section II. Obviously, (fc) = 2m. 
Fig. 4(a) shows the order parameter a vs m for differ- 
ent period T of the infection cycle. For a fixed period 
T, a transition in the synchronization can be observed as 
m increases. Moreover, the larger the period T, the less 
the critical value of mc, at which the transition occurs. 
We set T — 10 to analyze the critical scaling by using 
standard finite-size analysis mentioned above. Fig. 4(b) 
displays that when (3/v = 0.16 the curves with different 
sizes N give a unique crossing point at m c ~ 2.8 (where 
(fc) c w 5.6). In Fig. 4(c), the fit gives (1 - B)/V ~ 0.35. 
Hence B ~ 0.31 and v ~ 1.95. 

V. Conclusion. — To summarize, we have investigated 
the influence of the strength of community structure (Q) 
on global and local synchronization induced by the SIRS 
epidemic dynamics. The numerical results have shown 
that small Q induces better global synchronization and a 
phase transition occurs at Q c ks 0.835 estimated by using 
finite size analysis; While for the local synchronization 
there exists a minimal value of order parameter a around 
Q w 0.75. This result is in accordance with Ref. 23] in 
which a modified simulated annealing algorithm is ap- 
plied to optimize the synchronizability and well-defined 
communities do not exit in the emerging networks. That 
implies the networks with small Q are of strong synchro- 
nizability. It is also worth mentioning that, as in syn- 
chronization process well-defined communities of nodes 
emerge in different time scales, Arenas et al have used 
the synchronization to reveal the community structure 

M- 

Moreover, we have studied the synchronization order 
parameter vs (k) on scale-free networks with Q = 0. The 
simulation results demonstrate that, for a fixed period 
T, a transition in the synchronization can be observed as 
(k) increases. The larger period T corresponds to smaller 
critical value of transition point (k) c . 

We acknowledge the support from the National Nat- 
ural Science Foundation of China under Grants No. 
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